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ABSTRACT 


A geometrically similar family of three supercavitating 
hydrofoils was tested in the M.I.T. variable pressure water 
tunnel. The half-span foils were of elliptical planform; 
the ratio of foil half-span to tunnel height was 1/4 for the 
small foil, 1/2 for the medium foil, and 3/4 for the 
largest foil. The tunnel was of square cross section. Lift, 
drag, moment, tunnel speed, ambient pressure, and cavity 
pressure were measured for attack angles from 8 to 21 degrees 
and a variety of ambient pressure settings; cavity length 
measurements were obtained from photographs. Results were 
compared with theoretical results obtained by Leehey, and 
also with a more detailed numerical lifting surface procedure 
developed by Jiang and Leehey. For the small and medium foils, 
it was sufficient to correct only for the effect on downwash 
of the images of the trailing vortices. The large foil data, 
however, required further correction; upon application of the 
same corrections which were applied to the data for the two 
smaller foils, the force and moment data for the large foil 
plotted slightly lower than did the data for the two smaller 
foils, while the cavity length data for the large foil 
indicated cavity lengths significantly larger than for either 
of the theoretical predictions or the cavity length data for 
the two smaller foils. Through the application of existing 
two-dimensional corrections, the force data for the large 
foil were brought into close agreement with the force data 
for the two smaller foils, but no suitable correction factors 
exist for the cavity length data. 
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NOMENCLATURE 


foil leading edge bevel length; see Fig. 2 
aspect ratio = 2s*/s 

foil mean chord, measured at centroid position 
foil root chord 


tunnel correction factor for influence of images 
of trailing vortex system 


drag coefficient = totai drag = Cy + 5E,: 
(1/2) pu2s_- uncor i 


uncorrected drag coefficient, as calculated from 


the raw experimental data = measured crag 
(1/2) pu“s 


drag coefficient, corrected for blockage; see 
equation (1) 


change in induced drag due to effect of images 
of trailing vortex system 


measured lift 


lift coefficient = (1/2) 002s 


two-dimensional lift coefficient; see equation 
(II-5) 


moment coefficient about the midchord 


point = moment 
(1/2) pU2Se 
tunnel height = 51 cm 


cavity length as a function of spanwise position; 
measured from leading edge 


cavity length at centroid position; measured 
from midchord 


minimum pressure along tunnel wall in vicinity 
of cavity 


q(x,2z) 


Measured cavity pressure 


upstream tunnel static pressure at the spanwise 
position of the foil centroid 


vapor pressure of water 

tunnel static pressure controller setting 
static pressure Signal to controller from 
upstream, downstream tunnel static pressure 
taps; see Fig. 3 


strength of discrete source located in the 
(i,j)-th element 


source strength at (x,z) 

foil half-span 

element length for (i,j)-th element 
planform area of half-span model 

tunnel test section cross-sectional area 
thickness of hydrofoil 


perturbation velocities in the x,y-directions 


water velocity one meter upstream of dynamometer 
shaft axis 


Maximum water velocity (at point of minimum 
pressure, Pp) in vicinity of cavity 


chordwise coordinate of leading edge as a 
function of spanwise position 


distance of cavity pressure tap from midchord 
foil coordinates; see Fig. 2 


distance of centroid of half-span foil model 
from midspan 


distance of cavity pressure tap from midspan 


geometric angle of attack with reference to 
tunnel side wall; when used in Appendices II 
and III, a indicates angle of attack in a free 
stream 


"true" angle of attack, as corrected for effects 
of images of trailing vortices = a + bo. 


change in angle of attack due to effect of 
images of trailing vortices 


arctan (2a/o,) 


strength of discrete vortex located in (i,j)-th 
element 


vortex strength at (x,z) 


ratio of foil frontal width to tunnel width; 
see equation (3) 


ratio of foil frontal area to tunnel cross- 
sectional area (without foil); see equation (4) 


dummy variables in the x,z-directions 
density of water 


cavitation number based on measured cavity. 


pressure = (1/2) U2 


cavitation number, based on vapor 


pressure = —Ea-Ei- 
(1/2) PU 


cavitation number, corrected for blockage; see 
equations (2) and (5) 


cavitation number, based on Pps Por and V, 


the maximum water velocity = —Pb-Pc_ 
(1/2) ove 
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INTRODUCTION 


Corrections for water tunnel wall effects during the 
testing of subcavitating hydrofoils are essentially the 
same corrections developed fifty years ago for use in wind 
tunnel testing of airfoils; references [1] and [2] are 
standard texts in this area. Only within the last twenty 
years have wall effects been considered in the case of 
supercavitating hydrofoils, with the majority of theoretical 
and experimental work concentrated on two-dimensional hydro- 
foils. The problem of blockage can be significant; the 
shape of the foil and cavity combination must be known in 
order to apply a blockage correction of the type utilized in 
subcavitating flow. Unfortunately, this proves to be an area 
of great difficulty in the various theories. Corrections to 
drag coefficient and cavitation number for the two-dimensional 


pure-drag case are given in references [3] through [7]; a 


nonlinear theoretical model of the two-dimensional lifting 
problem is formulated in reference [6]. In reference [8], 
Baker converts this model into a computer program which 
numerically calculates the wall effect on two-dimensional 
hydrofoils of arbitrary shape. It was hoped by Baker that 
these results wouid prove adequate for application to flows 
over high aspect ratio hydrofoils. Prior to the current 


research, a systematic series of experiments which exhibit 
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wall effects on a geometrically similar family of hydrofoils 
had not been carried out in either the two- or three- 
dimensional case. The purpose of the present research is 

to make available data on a geometrically similar family of 
three supercavitating hydrofoils of finite span, and to apply 
existing two-dimensional corrections in an attempt to bring 
the data into agreement with three-dimensional, unbounded 


flow theory. 


nt Fe 


EXPERIMENTAL APPARATUS AND PROCEDURE 


Apparatus 


The experiments were conducted in the M.I.T. recircula- 
ting, variable pressure water tunnel (Fig. 1). The test 
section is 51 cm square in cross-section, and 147 cm long. 
Test section velocity was continuously variable between 0 
and 9.1 m/sec; a manometer indicated the difference in 
pressure across a contraction in the tunnel one meter up- 
stream of the dynamometer axis. This manometer had been 
previously calibrated using a pitot tube in the test section. 
The assumption was made that the blockage caused by the foil 
and cavity combination did not affect the manometer 
appreciably. Test section static pressure was variable 
between 60 mm Hg absolute and atmospheric pressure; static 
pressure readings were taken using a mercury manometer 
connected to either or both of two taps on the centerline 
of a side wall of the test section. These taps were 
located 62 cm upstream and downstream of the hydrofoil 
dynamometer, which was located at the center of the test 
section. 

Foil dimensions are given in Fig. 2. The three half- 
span, aspect ratio five foils had sharp leading edges. 
Geometrically similar elliptical planforms were chosen even 
though foils of this type were more difficult to fabricate 


than, for example, foilsof rectangular planform; elliptical 


at 


planform foils exhibit less tip loading effect than do 
rectangular planform foils, thereby rendering them more 
amenable to analysis. Each foil was drilled with a hole 
terminating in a pressure tap on the cavity side of the 
foil. The medium-sized foil was the aspect ratio five foil 
utilized by Leehey and Stellinger [9], except that the 
cavity pressure tap was relocated to a position nearer the 
tip. The composition of the medium foil was of type 304 
Stainless steel, while that of the small and large foils 
was of type 614 aluminum bronze. Each foil was welded to 
its own 10.16 cm diameter mounting plate; the mounting 
plates, in turn, were bolted to the dynamometer so that 
the small, medium, and large foils protruded 1/4, 1/2, and 
3/4 of the way, respectively, downward from the top window 
into the test section. 


The dynamometer was capable of measuring three 


moment and three force components, but only the moment 

about midchord and the forces tangential to and perpendicular 
to the span and chord were measured. Zero geometrical 

angle of attack was set by aligning the flat (wetted) side 

of each foil with the test section side wall. Upon 
alignment, a sighting telescope, which was attached to the 
dynamometer, was aligned with a scale on the laboratory 

wall; this scale, graduated in increments of 0.05 degree, 


was utilized in setting the geometrical angle of attack. 
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| Procedure 

Before test runs were begun on any one day, the tunnel 
was operated at medium velocity for twenty minutes at a 
pressure of 300 mm Hg absolute; this removed almost all of 
the air trapped in remote locations of the tunnel. Prior 
to each test run, the bubble chamber shown in Fig. 3 was 
utilized to collect and draw off any gases or air which may 
have been trapped in the static pressure sensing line; once 
Cleared of gases, the bubble chamber was merely a reservoir 
of common pressure. On the other hand, small amounts of air 
were bled intermittently into the cavity pressure sensing 
line during testing in order to keep it free of moisture. 

Each foil was tested at a variety of attack angles 
from 8 to 21 degrees. The lower limit of 8 degrees was 
chosen so that the cavities would not terminate on the foil. 
For all test runs, the test section velocity was set at 


8.5 m/sec, whereupon static pressure was reduced until a 


cavity covered the entire planform. Readings were taken 

of lift, drag, moment, static pressure (with pressure taps 
U and D open; see Fig. 3), cavity pressure, and velocity. 

A Polaroid photograph was taken of the foil and cavity from 
the wetted side of the foil, utilizing an exposure time 

of 1/40 second with bottom flood lighting. Pressure tap D 
was then closed, and the same readingstaken; reopening 
pressure tap D and closing pressure tap U, the same series 


of readings was repeated once more. (Photographs were not 
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taken for these last two pressure tap combinations since 
hardly any variation in cavity shape or length was noted 
when compared to the situation where both pressure taps 
were open.) This completed the measurement of three data 
points at one controller setting. Both pressure taps 

were then opened, and the static pressure was reduced by 
20 to 30 mm Hg, whereupon the above procedure was repeated. 
The static pressure was reduced in this fashion until an 
excessive number of recirculated cavitation bubbles resulted 
in either extremely poor conditions for photography or 
rapidly fluctuating instrument readings. Prior to and 
immediately following a run of this type, the room and 
water temperatures were recorded, as were the tares on the 
instrumentation. 

A computer program was utilized in the data reduction. 
Wetted surface frictional drag (for each foil and its 
mounting plate) and dynamometer shaft twist were taken into 
account, and test section static pressure readings were 
modified for static head difference to give the static 
pressure at the vertical position of the foil centroid. 
(Gravity effects were otherwise neglected.) Cavitation 
number was computed using measured cavity pressure, and 
also by using vapor pressure at the ambient water temper- 
ature (computed using equation (2) on page 151 of 


reference [10]). Lift, drag, and moment coefficients 
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were also calculated. Standard wind tunnel corrections 
(as described in Appendix I) were applied to the data to 
account for the effect of the images of the trailing 
vortices; this resulted in an increase in drag coefficient 
and also in an increase in effective angle of attack. 
Utilizing the photographs, cavity lengths were measured 


from the midchord at the centroid position. 


: 
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RESULTS 


Experimental results were compared with predictions 
from Leehey's linearized theory [11], which utilizes the 
method of matched asymptotic expansions; the theory is valid 
to first order in angle of attack, and to second order in 
the reciprocal of aspect ratio. Comparisons were also 
made with a more detailed linear lifting surface calculation 
developed by Jiang and Leehey [12]. This numerical theory 
is expected to be more accurate at the lower values of 
0 /Ons where the cavity length is the same size as the foil 
span or longer; it should also be more accurate, in general, 
for prediction of moment coefficients. Synopses of the 
theories developed in references [11] and [12] are given 
in Appendices II and III, respectively; the experimental 
data are given in Appendix IV, both in raw form and with 
standard wind tunnel corrections (Appendix I) applied. 

The wall boundary layer momentum thickness was 
approximately 1.9 mm, as inferred from measurements by 
Stellinger under similar conditions [9]. Since less than 
two percent of the wetted surface area of the small foil 
was within 1.9 mm of the upper wall, the upper wall 
boundary layer was assumed to have negligible effect on the 
force and moment coefficients of all three of the foils 


tested. 
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It should be noted that the only correction for wall 
effect which has been applied to the data in Figs. 4 
through 30 is the correction for the effects of the images 
of the trailing vortices (Appendix I). Also worthy of note 
is the fact that the supercavitating condition was never 
achieved for the small foil at an angle of attack of 8 
degrees. This fact bears upon the lack of scaling evident 
in the force and moment data of Figs. 4, 12, and 20 for this 
foil. This matter will be considered again in the sub- 


sequent discussion section. 


Lift (Figs. 4 through 11) 


There was very good agreement, overall, between 
theoretical predictions and the experimental data. For 
small values of the similarity parameter F/O (i.e., for 
long cavities), there was much better agreement with the 
numerical theory than with the asymptotic theory; in this 
case, the foil and cavity combination no longer has a 
large aspect ratio. For the higher values of 0 /One 
Leehey's theory appears to underpredict, especially at the 
higher angles of attack. This is in apparent contradiction 
with Fig. 4 of reference [9]; it should, however, be noted 
that the data of reference [9] were based on cavitation 


numbers calculated with vapor pressure rather than with 


measured cavity pressure. It can be seen from Figs. 32 
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through 35 of the present work that, when based on measured 
cavity pressure, the cavitation number will be smaller, so 
that the above contradiction is only apparent. The small 
and medium foil data plotted almost on top of one another, 
whereas the large foil data plotted lower than did the small 
and medium foil data, especially at attack angles greater 


than 11 degrees. 


Drag (Figs. 12 through 19) 

When the effects of streamwise foil or flow curvature 
are negligible, linear theory predicts Cy = OnCy so that 
the drag data are plotted as Cp/an vs. T/A Ch/an plotted 
somewhat higher than the theoretical predictions except 
for very long or for very short cavities. As with the lift 
data, the small and medium foil data plotted together, with 


the large foil data plotting somewhat below these, the 


effect becoming more pronounced at the higher angles of 


attack. 


Moment (Figs. 20 through 27) 


The moment coefficient was taken about the midchord 
consistent with the right-hand rule. As pointed out by 
Leehey and Stellinger [9], Leehey's matched asymptotic 
expansion theory neglects lifting surface corrections 
(third order in the reciprocal of aspect ratio) for the 


moment coefficient. It is not surprising, therefore, that 
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the current experimental data agreed best with the more 
detailed lifting surface theory of Jiang and Leehey [12]. 
The load cell utilized for moment measurements on the small 
foil was of inadequate sensitivity, and consequently, much 
of the small foil data at the smaller attack angles may be 
less accurate than the data for the medium and large foils; 
comparisons must therefore be made between the large and the 
medium foil data. Once again, the large foil data plotted 
below that of the medium foil, especially at the larger 


attack angles. 


Cavity Length (Figs. 28 through 31) 


As mentioned previously, cavity lengths were measured 
from midchord at the centroid position rather than from the 
leading edge, since this convention was utilized by Leehey 
Lil]. Cavity length (L) was nondimensionalized on the foil 
mean chord (c). The cavity length data for the small and 
medium foils plotted on top of one another (Figs. 28 and 29); 
in both cases, the cavity lengths were slightly less than 
those predicted by theory, but the overall agreement with 
theory was good. In contrast, there was a marked deviation 
in the large foil cavity length (Fig. 30) and shape (Fig. 31), 
especially at the root, where local blockage was the worst. 
The tip vortices seen in the photographs had no evident 


effect upon the force or moment coefficients. Neither 
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Leehey's matched asymptotic expansion theory nor the numer- 
ical lifting surface theory developed by Jiang and Leehey 
accounts for the roll-up effects which produced these 
isolated trailing vortices in the wake. 

Representative photographs of the large, medium, and 
small foils at the same angle of attack and at approximately 
the same cavitation number are shown in Fig. 31. The cavity 
shapes shown were averaged over the exposure time of 1/40 
second; visual observations showed that the instantaneous 
cavity shape was highly unsteady. For the small and medium 
foils, the general outline of the cavity was elliptical, 
as predicted by theory; the cavity length was undoubtedly 
reduced toward the tip due to the hydrostatic pressure 


gradient existing in the tunnel. 


Cavitation Number (Figs. 32 through 35) 


Below a cavitation number of approximately 0.4, there 
was almost no difference between cavitation number calculated 
with measured cavity pressure or with vapor pressure. It 
should be noted that in this region there exists a large 
variation in cavity length from very short to very long. 

This indicates that for even short cavities, where ram 
effects of the re-entrant jet from the cavity end have 
heretofore been thought to have a pronounced effect on 
measured cavity pressure, the actual cavity pressure was 


approximately equal to the vapor pressure of the water. 
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Above a cavitation number of 0.4, the disparity between 
a and o, increased with increasing cavitation number; the 
data in this region, however, were only for large angles of 
attack and small cavity lengths, indicating that there may 
indeed have been some effect caused by the impinging of the. 
re-entrant jet on the cavity pressure tap. 

After the series of basic experiments was completed, 
ram effects were further investigated on the large foil 
by Jiang and Leehey, this time utilizing a total head tube 
for cavity pressure measurement. The L-shaped total head 
tube protruded downwards into the cavity from the upper 
tunnel wall so that it was parallel to the foil wetted 
surface and pointed towards the leading edge and away from 
the impinging re-entrant jet. Cavity pressure readings 
were taken with both the standard foil surface pressure tap 
previously utilized and the total head tube, and vapor 
pressure was calculated as was done previously; results are 
shown in Fig. 35. For the larger attack angles and shorter 
cavities (higher cavitation numbers), use of the total head 
tube significantly reduced re-entrant jet effects; the 
result was that cavitation numbers calculated with total head 
tube data were as much as 5 to 8 percent higher than 
cavitation numbers calculated with data obtained from the 
previously used foil surface tap. Although a similar check 
was not made on the small and medium foils, it can be 


assumed that the same general trend would occur. It should 
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be noted, however, that data points in Figs. 32 and 33 


exhibit a more pronounced downward "hook" at higher 


cavitation numbers than do the large foil data (Figs. 34 


and 35), a fact which accounts for the noticeable "hooks" 
in the lift and moment data at the higher angles of 

attack (see, for example, Figs. 10, 1l, 26, and 27). 

If the total head tube had been utilized throughout the 
entire series of experiments, these deviations from theory 
at high attack angles and short cavities would probably 


have been less. 


Influence of Static Pressure Controller 

As stated earlier, readings were taken for the various 
combinations of the two static pressure taps; this was 
done in order to determine the effect of blockage on 
pressure and velocity at the two different tap positions. 
It was anticipated that this variation in static pressure 
tap combination would have no effect on the forces and 
moments experienced by the foil. Due to the location within 
the system of the static pressure controller (Fig. 3) 
(which was inadvertently overlooked), the forces and moments 
varied significantly when the combination of static pressure 
taps was changed. The explanation for this is as follows. 
Before data were taken, both static pressure taps were opened 
and the controller setting reduced to 5; the actual static 


pressure in the tunnel was decreased by the control system 


eal 
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until the combined signal from both taps was p. When the 
instrument readings settled out, data were taken. The 
downstream static pressure tap was then closed, so that the 
only signal to the controller was Py: which was greater than 
the controller setting, p. The controller therefore lowered 
the actual static pressure in the tunnel until Py equaled Ds 
Consequently, the effective cavitation number (as seen by the 
foil) decreased, causing the cavity to grow. This reduced 
the effective camber of the foil and cavity combination; 

this loss of effective camber resulted in a reduction of the 
force and moment coefficients from their previous values. 

On the other hand, the measured cavitation number increased 
slightly due to a small decrease in measured cavity pressure; 
the presence of the controller prevented the measured static 
pressure from changing significantly. The opposite effect 
was observed with the upstream tap open and the downstream 
tap closed ‘i.e., an increase in force and moment coefficients, 


an increase in effective cavitation number, and a slight 


decrease in measured cavitation number when compared to the 
case when both pressure taps were open); there were only a 
few of these data points taken. 
Although the desired pressure readings were not obtained 
in the current experiment (i.e., the actual tunnel static 
pressure was cifferent for each of the tap combinations), 


the results do show the significance cf the blockage. When 


aga 


comparing the results obtained with the different tap 
combinations for the small and medium foils, there was very 
little change in the plotted data points; the increased 
blockage in the case of the large foil resulted in signi- 
ficantly different force and moment data. For the sake 
of clarity, the data for the different tap combinations are 
shown only for the large foil at an angle of attack of 
| 14 degrees. Unless otherwise specified, further remarks 
concerning foil data will be for the data obtained with both 
static pressure taps open. 

The general trend of the data agrees with the experi- 
mental results obtained by Kermeen (See page 37 of reference 
[13]); he observed that the longer the cavity in super- 
cavitating flow, the smaller the force coefficients. The 
same conclusion follows from the theoretical results of 
Leehey, as may be seen from Fig. 4 of reference [9]. 

At a later time, the system indicated by the dashed 
lines in Fig. 3 will be utilized to obtain data on the large 
foil; the present system would be operated with only tap U 
open to provide an indication of static pressure at "infinity" 


to the controller. This configuration will allow pressure 


measurements to be made without affecting the input signal 


to the static pressure controller. 


 —— 
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DISCUSSION 


For the small and medium foils, it appears to be 
necessary to correct only for the effects of the images of 
the trailing vortices in accordance with standard wind 
tunnel procedure. These corrections to drag coefficient 
and effective angle of attack brought the data into good 
agreement with theory; the lift, moment, and cavity length 
data also agreed well with theoretical predictions. For 
the most part, the plots for the small and oan foils 
were very close to each other. 

Upon application of the same corrections to the large 
foil data, good agreement with theory was seen for the force 
and moment coefficients; they were, however, generally 
lower than the data for the smaller foils, especially at 
higher attack angles. The plot of cavity length for the 
large foil, however, showed substantially longer cavities 
than predicted by either of the two theoretical models used 
for comparison. The most pronounced deviations in cavity 
length occurred for the long cavity (lower F/O) data. 

When compared with the data obtained from the two 
smaller foils, the general trend of the large foil data is 
in agreement with Baker's predictions for wall effects on 
two-dimensional supercavitating hydrofoils [8]. Baker makes 


the point that the wall effect is negligible for tunnel 


LL 


ly 


height-to-foil chord ratios (H/c) greater than 10 for 
cambered foils. Although he makes no such generalizations 
concerning supercavitating flat plates, (i.e., the case 
applicable to the current research), it is seen from 

Figs. 4 through 27 that the effect of blockage on the force 
and moment coefficients is not exceedingly large. The 
following line of reasoning explains why the force and 
moment data for the large foil were somewhat lower than for 
the two smaller foils. For the three foils at the same 
cavitation number and at the same angle of attack, the 
magnitude of the velocity on the cavity wall of all three 
foils is 9,/2 (from linearized theory); consequently, all 
three foils have identical pressure coefficients on the 


cavity. The largest foil experiences significantly more 


blockage than does either of the smaller foils, so that the 
velocity on the wetted surface of the large foil is much 
higher than for either of the smaller foils; the pressure 
coefficient on the wetted surface of the largest foil is 
therefore the least of all three foils, resulting in smaller 
lift, drag, and moment coefficients than for either of the 
smaller foils. It is interesting to note that this overall 
effect is contrary to blockage effects experienced in 
subcavitating water tunnel or wind tunnel testing, where 


force and moment coefficients are increased. 
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Baker also questions the sharp increase in wall effect 
predicted .; Wu, Whitney, and Lin [6] for thin symmetric 
wedges or for small attack angles in the lifting flow case. 
The results of the present work do not show a marked 
increase in wall effect on force and moment coefficients at 

: the lower angles of attack, in agreement with Baker. 

Conversely, Baker's results indicate that the wall 
effect on the cavity plus wake can be tremendous. In the 
current research, this is, in fact, the most significant 
departure of the large foil data from the theoretical pre- 
dictions and from the data of the two smaller foils (Figs. 28 
through 30). 

Several attempts were made to bring the large foil data 
into closer agreement with the small and medium foil data by 


using existing two-dimensional blockage corrections. Wu, 


Whitney, and Brennen [4] derive wall effect corrections 
for the two-dimensional pure-drag (wedge) case, employing 
the open-wake and Riabouchinsky models. Wu, Whitney, and : 
Lin [6] consider wall effects in two-dimensional lifting | 
cavity flows, a situation more applicable to the current 

research; Baker [8] utilizes a computer program to translate 

this into information which could be readily applied to 


correct experimental data. 
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The pure-drag, open-wake model employs corrections to 
both drag coefficient and cavitation number. The drag 
coefficient corrected for blockage, C,,', can be expressed as 


l+o' 


2 
Cy + O(A°), (1) 


sili 


where o', the cavitation number corrected for blockage, is 


given by 


l+c | 2 
| cyA + O(dA), 


% 


where o_ is the measured cavitation number, C 


é is the drag 


D 
coefficient as defined in the nomenclature, and i is the 
ratio of foil frontal width to tunnel width, which may be 
expressed as 

(mean chord) x sin(a) 


=:—_—_—_ (3) 
tunnel width 


For the large foil with \ defined in this fashion, \ = 0.30 
gin (a). Application of the above corrections to the large 
foil data shifted the data points down and to the left on 

the plots of C)/as vs. o(/a, (Figs. 12 through 19) but did 
not bring the data into any closer agreement with the small 


and medium foil data. When the cavity length data for the 
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large foil were replotted using the corrected cavitation 
numbers, the data points in Fig. 30 were moved downward; 
this downward movement was not sufficient, however, to bring 
the large foil cavity lengths into agreement with either 
theory or the small and medium foil data. 

Actually, the straightforward application of such a 
two-dimensional correction to the case of a finite-span 
wing should overcorrect for blockage, making A much larger 
than it should be. Perhaps a better representation of A 
for the finite-span case would be 


foil frontal area 


a Sa a 
3D tunnel cross-sectional area (without foil) 


(4) 
(foil planform area) x sin (a) 


ee oo Ed 
= a 


tunnel cross-sectional area (without foil) 


Recalculating the corrections using i in place of id in 


3D 
equaticms (1) and (2) would make the corrections even smaller 
than previously calculated. Since the previous corrections 
proved to be inadequate, there is no reason to believe that 
the modified corrections should be any better. 

To correct measured cavitation number and measured 
drag coefficient utilizing the Riabouchinsky model, it is 


necessary to obtain Pp’ the minimum pressure reading along 


the wall in the region of the cavity. Equation (1) is again 


a le 


used to correct the drag coefficient, but now the corrected 


cavitation number is calculated as 


o'= (2/3)0,+ (1/3) 0" 


where o'' is the cavitation number based on Por Pp: and 


maximum velocity, V, and where O, is the cavitation number 


based on Por Pos and upstream velocity, U. It had orginally 


been hoped that a reasonably close measurement of the minimum 


pressure could be obtained by using only the downstream tap 


(tap D in Fig. 3); as was shown earlier, the interaction of 


the static pressure controller precluded making this measure- 


ment accurately. In an attempt to apply the correction based 


on the Riabouchinsky model, the gross assumption was made 


that the maximum frontal area of the cavity was approximately 


the same as the frontal area of the foil; utilizing a 


continuity argument and the steady-flow Bernoulli equation, 


the pressure minimum was estimated, thus allowing application 


of the wall effect correction based on the Riabouchinsky 


model. Unfortunately, this procedure gave results which 


were comparable to those obtained using the open-wake model. 


The failure of both models in the current application 


should not be surprising since they were designed for two- 


dimensional, pure-drag cavity flows, and not for three- 


dimensional, lifting cavity flows. 
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Whereas the pure-drag wall effect corrections changed 
both drag coefficient and cavitation number, the results 
of Baker's computerization of the two-dimensional Wu, Whitney, 
and Lin [6] lifting case wall effect model are presented 
in terms of a change to the force coefficients (i.e., 


drag and lift coefficients in this case), referenced to the 


uncorrected measured cavitation number, TQ: Accordingly, 


Baker's corrections to the lift and drag coefficients are 
expressed as a percentage difference from the free stream 
condition. He presents data for only one ratio of tunnel 
height-to-arclength (i.e., chord length), H/c = 4. 
Calculating this ratio using the mean chord of the large foil 
gives H/c = 3.33; recalculating H/c as the ratio of tunnel 
cross-sectional area to foil area, H/c = 4.48 is obtained. 
It was concluded, therefore, that the H/c = 4 corrections 
given by Baker should be suitable for correction of the large 
foil data. Upon application of these additive corrections 
to both lift and drag data, it was found to bring the large 
foil data into close agreement with the small and medium 
foil data. 

Because an infinite wake model is utilized in 
reference [8], Baker's results give no correction for the 
extremely large cavity lengths recorded for the large foil. 


His results do show, however, that as blockage increases, 


the nondimensional cavity thickness increases. This was 
borne out nicely when the cavitation characteristics of the 
three foils utilized in the current research were compared 

at an angle of attack of 8 degrees. On the small foil, 

which produced the least blockage, supercavitation was never 
totally achieved, no matter how low the static pressure was 
brought. For the medium and large foils, which produced more 
blockage than did the small foil, cavity termination was 


achieved behind the foil over a range of static pressures. 
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CONCLUSIONS 


As a general rule, the agreement of the experimental 
data with unbounded flow theory was good for the small and 


medium foils, whose ratios of half-span to tunnel width 


= aah RN 


were 1/4 and 1/2, respectively, and whose ratios of tunnel 
height to mean chord were approximately 10 and 5, respectively. 
It appeared that in order to bring the small foil and medium 
foil data into agreement with theory, it was necessary to 
apply only standard wind tunnel procedures in correcting 
for the effects of the images of the trailing vortices. 

of the three geometrically similar aspect ratio five foils, 
the largest foil, which protruded 3/4 of the way into the 
tunnel, and which had a ratio of tunnel height to mean chord 
of approximately 3.33, had force and moment data which were 
close to the small and medium foil data, although slightly 
lower; the cavity length data for the large foil differed 
greatly from the small and medium foil data, showing much 
longer cavities for the same ratio of cavitation number to 
angle of attack. It is evident that blockage corrections 
are required in addition to the standard wind tunnel down- 
wash corrections. 

Application and modification of blockage corrections 


based on various two-dimensional pure-drag cavity flow models 


failed to bring the large foil data into agreement with the 
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data for the two smaller foils. It was found, however, that 
corrections based on a two-dimensional lifting cavity flow 
model brought the large foil force coefficient data into 
much closer agreement with the small and medium foil 
data. While these force coefficient corrections were found 
to give an "engineering order of accuracy," no corrections 
have been found which give adequate corrections to the large 
foil cavity length data; it is clearly evident that further 
analytical work is necessary in this area. 

It is recommended that for future experiments, a static 
pressure measuring system similar to that depicted in Fig. 3 


be utilized. 


10. 
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(b) MEDIUM FOIL: @=12°, Of: 3473, Le+155 
¢ 


(C) LARGE Faw. +12", Of: 3268, L/c# 2! 


FIG. 31: COMPARISON OF CAVITY LENGTHS 
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APPENDIX I 
Standard Wind Tunnel Corrections for Effects of 
mages of Trailing Vortices 

The experimental data shown in Figs. 4 through 30 are 
corrected for the apparent upwash due to the interference 
of the tunnel walls, in accordance with standard wind 
tunnel practice. The boundary condition of zero normal 
velocity at the tunnel walls is satisfied in the analytical 
treatment of the problem by introducing a doubly infinite 
system of image vortices; the interference due to the 


presence of the tunnel walls is then calculated to be the 


interference caused by the series of vortex images on the 
flow at the model. The net result is an additive correction 
to the induced drag and the angle of attack (to be perfectly 
correct, there is also a slight negative correction to the 
measured lift, but this is usually neglected). Upon 
application of the correction factor, the flow about the 
model is transformed to the flow which would be seen if the 
foil were in an infinite stream. 

Since the hydrofoils tested in the square tunnel were 
only of half-span, a full-span model would require a duplex 
tunnel (that is, one whose ratio of height to width is 1/2). 
The ratio of model span to tunnel width for the small, medium, 
and large foils is 1/4, 1/2, and 3/4, respectively. Entering 


Fig. 6.34 of Pope and Harper [1], &, the tunnel correction 


= Fhe 


factor for the influence of the images of the trailing 
vortex system was found to be 0.126 for the small, 0.104 for 
the medium, and 0.092 for the large foil. If the ratio of 
the foil area to tunnel cross-sectional area is S/S. then 


the change in angle of attack is 


da; = Ss- Cy: (radians) (I-1) 


The maximum value for Aa, encountered in the experiments was 
1.47 degrees. When Aa, is added to the measured geometrical 


angle of attack, a, the "true" angle of attack is obtained 
Op = a + Aa; - (I-2) 
(In addition to the above correction, the data reduction 


program accounted for the dynamometer shaft twist, which 


never exceeded 0.06 degree.) The change in induced drag is 


Ac 


ul 
(or) 


D; So L 


i] 
co 
eg 
a 


2 ae 


so that the total drag coefficient is 


G. =¢€ + AC (I-4) 
= uncor D; 


where Cy is the uncorrected drag coefficient calculated 
uncor 


from the raw experimental data. 


a15< 
APPENDIX II 


Leehey's Matched Asymptotic Solution for 
Supercavitating Hydrofoils of 
Large Aspect Ratio 


A synopsis of Leehey's theoretical results, as derived 
in references [9] and [11], is given. 

Leehey's linearized steady-flow theory, which utilizes 
the method of matched asymptotic expansions, is valid to 
first order in angle of attack, and to second order in the 
reciprocal of the aspect ratio. Gravitational effects were 
neglected in the derivation. The analysis shows that the 
lift, drag, and moment coefficients, and non-dimensional 
cavity length (measured from midchord) are given by the 


following relationships: 


TO. 2siny-1l 
Cc, = —_— Li - ———_}, (II~1) 


L  siny (1+siny) A(1+siny) 
na? 2siny-1 
Cy = ac, = ————_ {1 - ———— +}, (II-2) 
sin y (l+siny) A(l+siny) 
~4a 1-2siny+2sin2y 
Cy = PERE l1- ? (II-3) 
1(1+siny) A siny (1+siny) 
and 
8 a 2 
2 
L/c = sec’y - 1/2 - - Cy P (II-4) 
TaA Cc 2D 
where Cy , the two-dimensional lift coefficient may be 


expressed as 


=366 


TO. 
Ch = , (II-5) 
2D siny (l+siny) 
and where 
20 
Y = arctan |—|l, (II-6) 
% 
and 
Pa"Po 
cz-_—_—_— (II-7) 


© (1/2) pu2- 


Experiments performed by Leehey and Stellinger [9] showed 
that the above predictions for lift and drag coefficients 
and for cavity lengths were good; it was evident from the 
experimental results, however, that lifting surface corrections 


would be needed to improve the moment coefficient prediction. 


Py iy 


APPENDIX III 


Numerical Lifting Surface Theory for Supercavitating 
Hydrofoils of Finite Span 


A brief description of the numerical lifting surface 


theory developed by Jiang and Leehey [12] for supercavitating 


hydrofoils of finite span will be given. 


The steady flow theory is linearized for application 


to thin wings at small attack angles. Discrete sources and 


vortices are utilized in the representation of the physical 


model, and the coupled integral equations are reduced to a 


set of simultaneous algebraic equations. The cavity closure 


condition is applied stripwise on the cavity in an iterative 


process to obtain the desired cavitation number over the 


supercavitated planform. 


This numerical theory should prove to be more accurate 


than the asymptotic theory of Leehey [11]; an implicit 


assumption in the latter theory is that the cavity length is 


less than the span, so that for small d/o (i.e., long 


cavities), the numerical theory should indeed give better 


predictions. It should be noted, also, that for the 


T/A = 1.0 curves in Figs. 4, 5, 7, and 8 of reference [9], 


there is an ever-present difference between the experimental 


data and Leehey's asymptotic theoretical predictions for 


lift and drag coefficients. For long and short cavities, 


however, there is an even stronger disagreement between data 


«Fon 


and Leehey's predictions for moment coefficient (Figs. 9 
and 10 of reference [9]}. The numerical lifting surface 
theory was developed in an attempt to overcome these 
difficulties. 

The foil-cavity surface is collapsed to become a 
region in the (x-z) plane of the foil. The relationships 
for the jump in the tangential and normal components of 
the perturbation velocity upon crossing the foil surface 


and cavity projection surface are: 


u(x,z,*0) - u(x,z, 0) = -y(x,z), (III-1) 


v(x,2,'0) - v(x,z,0) = q(x,z), (III-2) 
where y and q are the vortex and source distributions, 
respectively. The boundary conditions are: 

dy (x,z, 0) on the foil 


Sitess Oe wetted surface (III~3) 


dx 
=! on the cavity = 
u= a /2 boundary ais 
& (2) 
o cavity closure 
| a(&,o)d& = 0 condition, (III-5) 
X_ (Zz) 


where %(z) is the cavity length, measured from the leading 


edge as a function of spanwise position, and X, (2) is the 


Gin 


chordwise coordinate of the leading edge as a function of 
spanwise position. 

The discrete vortex and source method is utilized in 
the formulation of a series of simultaneous equations. 
Upon subdividing the projected surface into small elements, 
a discrete bound vortex, a trailing vortex, and a source 
are located within each of these elements. The bound 
vortex is situated on the quarter-chord line of each element, 
and induced velocities are calculated for all elements at 
their midspan, three-quarter chord positions. The 
concentrated source is taken to be a constant distribution 
across each element at its three-quarter chord position. 
The boundary condition of fixed cavity pressure is imposed 
on each element at its midspan, quarter-chord position. 
Utilizing these boundary conditions, the following equations 


are obtained: 


pe ea p Md aa 
eee (III-6) 


ys se gM (III-7) 
LypsesgN 
= 0 for k > Ml 


LjpstegNy iiss) 


ijke are the v and u velocities at the (k,2) 


"ae a 
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~th control point induced by the unit strength vortex and 

source at element (i,j). For the (i,j)-th element, 484; is 
the element length, Vag 
the discrete source strength. Ml and N are the number 


the discrete vortex strength, and 


Wis 
of elements on the foil along the chord and span, respectively. 


M is the number of elements along the chord in the cavity. 
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APPENDIX IV 


: 
| Data Reduction Program and Data Printout 


Beginning on p. 82, a listing is given of the FORTRAN IV 
computer program utilized in the reduction of the experimental 
data. 

The program output is also given, beginning on p. 92. 

For each foil at each angle of attack, a complete printout 

of the experimental data is given in raw form, with corrections 
for instrument tares, in reduced form, in none inenataeat form 
without any corrections for wall effect, and finally in 
nondimensional form with corrections for the effects of the 


images of the trailing vortex system. 
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